Cooling classical particles with a microcanonical Szilard engine 
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The possibility of extraction of energy from a system in a cyclic process is discussed. We present 
an explicit example where a system, initially prepared in a microcanonical state, is able to perform 
such operation. The example is similar to the Szilard engine, but the microcanonical initial condition 
allows one to design a protocol where measurement is not necessary. 
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Extracting energy from a system faces limitations im- 
posed by Thermodynamics and Statistical Mechanics. 
Kelvin statement of the Second Law precludes the cyclic 
extraction of energy from a single thermal bath. Using 
Hamiltonian dynamics, Jarzynski [T] proved that this is 
the case for systems initially prepared in a Boltzmann 
state. Suppose a system described by a Hamiltonian 
H{q,p; Ai, . . . , A^), where (q^p) stands for a point in the 
phase space and Ai(t), . . . , Xn{t) are external parameters 
operated by an external agent. The system is isolated, 
except for the interaction with this external agent, and 
consequently the probability density p{q,p;t) obeys the 
time-dependent Liouville equation. Jarzynski pi proved 
that, along an arbitrary cyclic process A^(0) = A^(r), the 
average energy of the system necessarily increases, i.e., 
cyclic energy extraction is not possible. Campisi ^ fur- 
ther generalized this result proving that the extraction 
is impossible for any initial distribution p(g,p;0), such 
that the resulting probability density for the energy is a 
decreasing function. Both Jarzynski's and Campisi's re- 
sults can be considered as mechanical proofs of the Kelvin 
statement of the Second Law. 

On the other hand, the microcanonical distribution, 
where the total energy of the system is known with a 
given accuracy, does not fulfill Campisi's requirement. 
In fact, Sato [3 devised a particular example — a one- 
dimensional particle in a potential depending on one pa- 
rameter X{t) — where, starting in a microcanonical state, 
the energy of the particle decreases along a process with 
A(0) = A(r). Although this result is remarkable, the ex- 
traction of energy cannot be repeated, since after the first 
cycle the probability density is no longer microcanonical. 

In this Letter we present an explicit example where en- 
ergy can be extracted in a systematic way. It is a micro- 
canonical version of the Szilard engine [4 , consisting of a 
classical particle in a time-dependent potential. In con- 
trast with the original Szilard model, the microcanonical 
initial condition allows one to operate the engine with- 
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FIG. 1: (color online). Different steps in the microcanonical 
Szilard engine. At stage (b), we have plotted the particle in 
the left half but recall that, under the protocol described in 
the text, the particle may have enough energy to move all 
over the whole box. 



out measurements. Our example provides a new strategy 
for microscopic cooling of classical isolated particles, but 
also touches some fundamental issues regarding the me- 
chanical rationale of the Second Law. 

Our microcanonical version of the Szilard engine [4] 
consists of a single one-dimensional isolated particle of 
mass m in a potential modified by an external agent. 
The particle is initially confined in a box of length 2L 
and the potential is changed as shown in Fig. [l] In (a- 
b) , a barrier is raised at the center of the box up to height 
Vb at speed 7. The barrier potential is given by: 




\x\) 



\x\ < I 
otherwise, 



(1) 



for t G [0,^5/7]- In the second step (b-c), the right 
wall is moved from L to Lp = 2L -\- I with velocity a. 
Finally, along (c-d), the left wall is moved from —L to 
with velocity /3 and the potential recovers its initial 
configuration (but shifted to the right). 
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FIG. 2: (color online). Log-lin plot of the energy En versus 
the cycle number n. Parameters used are m = 1, L = 20, / = 
1.0, for different values of a, /3, 7. The maximum height of the 



barrier in cycle n is 



(n) 



^^^0, with K = [1/2 + //(6L)]' 



Now suppose that the barrier height Vb is equal to 
the initial energy of the particle Eq. When the barrier 
is raised along step (a-b), the particle gains some small 
amount of energy, and then the barrier is not high enough 
to confine the particle in one of the halves of the con- 
tainer. Along the expansion (b-c), the particle loses en- 
ergy in each collision with the right moving wall (in each 
collision the velocity changes from v to —v-\- 2a). After 
a certain number of collisions, its energy drops below Vb 
and the particle gets confined in the rightmost half of the 
container. As a consequence, the final compression (c-d) 
does not affect the particle much and the net effect of the 
process is a reduction of its energy. Since the particle will 
be always confined in the rightmost half of the box, no 
measurement is necessary along the step (b-c), unlike the 
original Szilard engine. This situation is possible only for 
certain values of Vb depending on the initial energy of the 
particle, although not on its position. We can therefore 
design a protocol that works for microcanonical initial 
conditions: fixed energy (velocity) and random position 
inside the box. Moreover, if the above steps are carried 
out quasi-statically, the energy change is deterministic. 
Then, by taking the appropriate values of Vb in each cy- 
cle, we can systematically extract energy at least down 
to a certain value which can be made arbitrarily small 
by reducing the speeds a, /3, and 7. Fig. |2] shows the 
energy after each cycle, following the described protocol, 
for different values of a, /3, 7. Data presented in Fig. |2] 
have been obtained by an event driving algorithm and 
solving exactly the Newton's equations for the motion of 
the particle. Along step (a-b), a third order polynomial 
must be solved for the motion along the barrier. 

For a full understanding of our microcanonical Szilard 
engine, it is convenient to study in detail the quasi-static 
limit, Q^,/3,7 0. Consider a point (xq^pq) in phase 



space, with position xq, momentum po? and energy £'0 = 
H{xo,po;0) =pg/(2m) + y(xo;0), evolving as {x{t),p{t)) 
with energy E{t) = H{x{t),p{t);t). The action 
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dx dp 



H{x,p;t)<E{t) 



(2) 



is an adiahatic invariant^ i.e., is constant for quasi-static 
changes of the Hamiltonian [5]|6]. In our case, the initial 
action (j)o = AL\/2mE(^ would in principle remain con- 
stant along the cycle. However, the invar iance has an 
exception: when an orbit changes abruptly due to segre- 
gation, confinement or sudden expansion, induced by the 
barrier, the action changes accordingly [6, 7 . Below we 
derive the conditions determining these changes. 

For a generic height V of the barrier and location of 
the left {Li) and right (L^) walls, the action of an orbit 
with energy E reads: 
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where or Li in the first line is taken depending on the 
side of the box where the orbit is confined. The discon- 
tinuity of the action (t){E) given by Eq. ([3| at £^ = V is 
a trademark of the breaking of its invariance. Suppose 
that the barrier is high enough to split, at some stage 
of step (a-b), the orbit of a particle with initial energy 
£^0- Right before the segregation, the energy £seg of the 
particle is equal to the height V of the barrier. Setting 
E = V in Eq. ([3| for £ > the invariance for the action 
(right before segregation) implies: 
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Segregation actually occurs if the maximum height Vb is 
greater than E^eg, i.e., if 



^0 < ( 1 - ^ ) Vb 



(5) 



Segregation decreases the action of individual orbits by 

A^seg = -<?^o/2 = -2L^^2mEo (6) 

Orbits lying in the right side of the box do not suffer 
any further collapse or expansion. On the other hand, 
those in the left hand undergo a sudden expansion when 
reaching an energy Vb along step (d-a). The increase of 
the action equals the action corresponding to a new lobe 
added to the orbit in the right side, which is given by the 
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right contribution in ^ a.t E = V = Vb {E < V) and 
Lj. = 2L + I, yielding: 



A0exp = 4v/2mVB [L 



(7) 



If there is no segregation, i.e., if the particle reaches 
stage (b) with an energy greater than Vb, then along 
step (b-c) the energy of the particle decreases, due to 
collisions with the right moving wall, and eventually can 
reach the critical value Vb confining the particle in the 
rightmost half of the container. In this case, the orbit 
of the particle suddenly loses the left lobe (see Fig. [s]). 
Before the confinement, the action is invariant. Setting 
E = V = Vb,Li = Lin Eq. ^ for E > V, the mvariance 
for the action, right before confinement, reads: 
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= 4Ly^2mEo. (8) 



We have written = Leonid the position of the right wall 
in the moment of confinement, which is at most Lp = 
2L -\- 1. Therefore, confinement occurs if: 
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Thus Eq. (jsj) and Eq. (|9|, constitute the condition for 
the confinement given by: 
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The change in the action due to confinement is the action 
corresponding to the left branch of the orbit disappearing 
after confinement, which is given by the left contribution 
in ^ dit E = V = Vb {E <V) and Li = L, yielding: 



'conf 
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Finally, if the initial energy is large, Eq > 
[3/2 + 1/{6L)]^ Vb, the orbit spans all over the box along 
the whole process and the action is invariant. Summa- 
rizing, the final action is 



■A(/>, 



(12) 



with = A(/)conf < if Eq. (10) is satisfied, Acf) = 



if the energy is above the upper limit in Eq. (10), and 
A^ = A^seg < or A^ = A^seg + A^exp > if it is 
below the lower limit. These two values are taken with 
probability 1/2, respectively. We can now calculate the 
final energy E^ since = + A^ = AL^J2mEf^ yielding 
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Fig. |3] shows the final energy as a function of the ini- 
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FIG. 3: (color online). Final energy as a func tion of the initial 
energy in the quasistatic limit, as given by ([l3|. Parameters 
used are m = 1, L = 20, / = 1.0. For small initial energies, 
particles are segregated to the left or right half of the box at 
step (a-b) with probability 1/2 corresponding, respectively, 
to the green and red curves depicted in the figure. The blue 
curve is the average final energy. 

tial energy. We see that the extraction of energy is only 
possible in the window defined by Eq. (10). However, 



since the energy change in the quasi-static limit is deter- 
ministic (if there is no segregation) we can tune Vb to 
verify Eq. (10) in subsequent cycles. The simplest choice 
is to set the height of the barrier at cycle n, equal 
to the energy of the particle when the cycle starts. In 
the case Vb = ^'o, the final energy after a cycle reads 
Ef = kEq with n = [1/2 + //(6L)]^ Therefore, an ap- 



propriate choice is Vj 



(n) 



^Eq, which is the protocol 



followed in Fig. [2] Notice that this choice is made a pri- 
ori and does not depend on the actual evolution of the 
system. Consequently, for finite velocities a, /3, 7 the pro- 
tocol fails, since the energy does not follow a determin- 
istic sequence as we depart from the quasi-static limit. 
However, the energy at which the protocol starts to fail 
can be made arbitrarily small by reducing the speed of 
the process, as Fig. |2] indicates. 

Our example prompts a crucial question: is the pos- 
sibility of extraction of energy a special feature of one 
dimensional systems or can be extended to systems 
with several degrees of freedom? A naive extension of 
the microcanonical Szilard engine to two or more in- 
dependent particles does not allow the systematic ex- 
traction of energy. Consider two one-dimensional par- 
ticles confined in the interval [— L,!/]. They are inde- 
pendent except for the initial preparation in a micro- 
canonical state with total energy Eq^ pEo{pi-,P2]^) = 
5 [Eq — p\/2m — P2/2rn) /(27rm). For simplicity, we will 
assume that the potential only affects particle 1. Then we 
can apply the above results using as initial distribution 
of energy the corresponding distribution of the energy of 
particle 1: 



p{E;0) = I dpidp25[E-j^]pEM,P2;0) 



Tr^/E{Eo - E) 



(14) 
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FIG. 4: (color online). Average change of energy after one 
(red) and an infinite number (black) of repetitions of the pro- 
tocol with fixed Vb and / = for two particles, as a function 
of the initial total energy Eq. Inset: probability distribution 
of the energy of the particle undergoing the protocol initially 
(filled blue), after one (red) and an infinite number of repeti- 
tions (black) for Vb = 1 and Eq = 2.5. 



Fig. [4] shows the change of energy for an infinitely nar- 
row barrier / = 0, as a function of Eq/Vb- The red curve 
corresponds to the change after one run of the protocol. 
One can extract energy in a single run for Eq > IAVb, 
but the probability distribution of the energy of particle 1 
changes in a rather uncontrolled way (see inset of Fig. |4j 
red curve), and a repetition of the protocol does not fur- 
ther decrease the energy. The black curve in Fig. |4] shows 
the change of energy after an infinite number of cycles 
(with the same Vb)- Adaptive protocols changing Vb in 
each cycle can slightly improve the energy decrease, but 
the systematic extraction of energy seems impossible in 
the case of two particles. With three or more particles, 
the probability distribution of the energy of particle 1 is 
a decreasing function and, following Campisi's theorem 
[2], the energy increases even after the first run of the 
protocol. 

However, the basic mechanism of our microcanonical 
Szilard engine is a "shuffling" of the phase space that 
could in principle work in systems with several degrees 
of freedom. Fig. [5] shows how the phase space changes 
after one run of the protocol for / = and a single 
one-dimensional particle. Points with energy below Vb 
(striped regions) are mapped to two regions of half size; 
those with energy between Vb and 9Vb/^ (blue regions) 
shift down, whereas the rest are not affected by the pro- 
tocol. Notice that this shuffling is fully compatible with 
Liouville theorem: the volume of any subset in the phase 
space is conserved. On the other hand, it needs the col- 
lapse or split of orbits. Otherwise, the adiabatic invari- 
ance of the volume enclosed by an energy shell implies 
that every point in the phase space goes back to its initial 
position after any quasistatic cycle. 

Up to our knowledge, there is no fundamental obstacle 



to reproduce this shuffling of the phase space in a system 
with many degrees of freedom. It is an open question to 
find an explicit example (probably more involved than a 
naive generalization of our protocol) or, on the contrary, 
to prove that the phenomenon described in this Letter is 
exclusive to systems with one degree of freedom. Either 
one or the other, the answer to this question touches 
a fundamental problem: the mechanical rationale of the 
Kelvin statement of the Second Law for systems prepared 
in the microcanonical state. Let us finally recall that the 
microcanonical ensemble, although of limited use in real 
applications, is essential for an objective formulation of 
Statistical Mechanics [8l[9]. 




FIG. 5: (color online). Evolution of regions in phase space 
after completing the protocol with a barrier height Vb and 
/ = 0, for a single one-dimensional particle. The striped 
region with energy lower than Vb, P ^ [—Pb,Pb], with 
Pb = V2mVB, splits into three regions leaving room to the 
blue region p G [ps, 3ps/2] to move down to p G [ps/2,ps]. 
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